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Abstract 

The dispersion properties of exciton polaritons in multiple-quantum- well based resonant photonic 
crystals are studied. In the case of structures with an elementary cell possessing a mirror symmetry 
with respect to its center, a powerful analytical method for deriving and analyzing dispersion laws 
of the respective normal modes is developed. The method is used to analyze band structure 
and dispersion properties of several types of resonant photonic crystals, which would not submit to 
analytical treatment by other approaches. These systems include multiple quantum well structures 
with an arbitrary periodic modulation of the dielectric function and structures with a complex 
elementary cell. Special attention was paid to determining conditions for superradiance (Bragg 
resonance) in these structures, and to the properties of the polariton stop band in the case when 
this condition is fulfilled (Bragg structures). The dependence of the band structure on the angle 
of propagation, the polarization of the wave, and the effects due to exciton homogeneous and 
inhomogeneous broadenings are considered, as well as dispersion properties of excitations in near- 
Bragg structures. 
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I. INTRODUCTION 



Optical properties of artificial structures with periodically modulated dielectric constant 
has been attracting a great deal of interest since pioneering papers Refs. Q and Q where 
such systems have been first discussed. Periodic modulation of the dielectric function signif- 
icantly modifies spectral properties of electromagnetic waves. Instead of a simple continuous 
spectrum with a linear dispersion law, the electromagnetic spectrum in such structures is 
characterized by the presence of allowed and forbidden bands similar to electronic band 
structure of crystals. For this reason, the new class of optical materials was dubbed pho- 
tonic crystals.— Changing the spatial distribution of the dielectric constant one can effectively 
control such fundamental properties as the group velocity of light, the rate of spontaneous 
emission, etc. This fact has important repercussion for both fundamental physics and for 
applications, where it opens up possibilities for new concepts of optical and optoelectronic 
devices. 

Modulation of the dielectric function, however, is not the only way to affect light prop- 
agation and its interaction with matter. It was shown in Ref. Ij that closely packed dipole 
active atoms can become coherently coupled by light, and this coupling significantly changes 
their emission properties resulting in the so called "super-radiance" effect. In the origi- 
nal Dicke's model^ the distance between the atoms was assumed to be much smaller than 
the wavelength of their emission, but similar effect can also arise if dipole active elements 
form a one- dimensional periodic lattice with the period coinciding with the emission half- 
wavelength (Bragg resonance). Such an arrangement is possible with a so called optical 
lattice of cold atoms^ and with their semiconductor analogs - Bragg multiple quantum well 
(MQW) structures.i^'2^ The latter are semiconductor heterostructures, in which very thin 
layers of a narrower band gap semiconductor (quantum well) are separated from each other 
by much thicker layers of a wider band gap semiconductor (barrier) in such a way that the 
period of the structure coincides with the half-wavelength of light emitted by excitons con- 
fined in a quantum well (QW). In these systems, excitons play the role of the dipole active 
excitations, which become radiatively coupled and demonstrate the superradiance effect. 
If, however, the size (the number of periods) in an optical lattice increases, the properties 
of the system chang e™ In particular, the dark modes form two branches of collective 
excitations, in which light is coupled with the material resonances, and which can be called 
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optical lattice polaritons. At the same time, the super-radiance mode develops into a stop 
band, which is a spectral gap between the two polariton branches, characterized by almost 
complete reflection of a normally incident radiation. The distinction between Bragg struc- 
tures and other arrangements, in which the period of the structure is not in the resonance 
with radiation emitted by the active elements, can, in this case, be described in terms of 
formation of a signiflcantly (by orders of magnitude) enhanced stop band. 

There is, however, an important difference between optical lattices of atoms and semi- 
conductor heterostructures. In the latter case, the interaction of light with active elements 
(excitons) is accompanied by multiple reflection of light from interfaces between wells and 
barrier layers caused by difference in their refractive indexes. Thus, in structures like 
MQWs optical lattice effects coexist with photonic crystal-like effects, which results in a 
number of new and interesting optical properties and opportunities for applications. Such 
structures, called resonant photonic crystals (RPC), have recently started attracting sig- 
nificant attention . ^^i^^i^^i^^i^^i'^'^'i'^^i'^^A general characteristic of this type of structures is the 
co-existence of dipole active material excitations (excitons, phonons, plasmons) described 
by optical susceptibilities of the resonant type and the periodic modulation of the back- 
ground dielectric constant. Two most fundamental problems in the focus of current research 
in this concerned with the effects of the interplay between resonances and spatial 

non- uniformity of the dielectric function on the band structure of these systems and their 
optical spectra. While these two questions are interconnected, they have to be recognized 
as two separate problems, one dealing with normal modes of closed (or infinite) systems and 
their dispersion laws and the other with the interaction of an internal radiation with a finite 
size structure. 

The main interest of studies in this field is a spectral region in the vicinity of the resonant 
frequency, where the interplay effects play the most important role. A theoretical description 
of RPC's in this spectral region is a challenging problem and its complete solution has not 
yet been obtained even for a simplest case of one-dimensional structures, such as MQWs. 
Of course, it is always possible to carry out numerical calculations of the optical spectra and 
the dispersion laws, which are particularly easy in one-dimensional case. However purely 
numerical approach does not provide a real understanding of physics of these structures. 
Therefore, it is very important to be able to carry out analytical analysis, which would not 
only provide a better understanding of physical processes taking place in these structures. 
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but would also be instrumental in designing structures with pre-determined properties, which 
is crucial for their possible applications. 

Different aspects of this problem have been, of course, considered in a number of previous 



publications. In particular, in Ref. 



201 the necessity of a modification of the Bragg resonance 



(super-radiance) condition for one-dimensional resonant photonic crystals compared to the 
case of optical lattice structures was shown. Later, this result was confirmed^i and an exact 
Bragg condition in such structures was found for the particular case of normal propagation 
of the electromagnetic wave. It was also noted in Ref. l2llthat. when the resonance condition 
is met, the spectral gap between the polariton branches becomes wider than in the case of a 
passive photonic crystal characterized by the same modulation of the dielectric function or of 
an optical lattice with the same strength of radiative coupling between the active elements. 
In the case of structures of higher (2 or 3) dimensions, no analytical results are available 



22, where it 



yet, but one can note recent numerical calculations presented in Refs. Q and 
was also found that the resonances may enhance polariton related band gaps compared to a 
purely passive photonic crystal with the same spatial distribution of the dielectric constant. 
It is important to emphasize that this enhancement of the band gap in RFC structures is not 
a trivial effect. In order to illustrate this point we can mention another, in a certain sense, 
opposite effect, which also results from the interplay between the resonances and periodic 
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inhomogeneities of the dielectric function. It was shown in Ref. 
frequency where the reflection coefficient of the incident light vanishes due to the destructive 
interference between the two channels of interaction between light and the structure. 

These examples show the richness of the optical properties of even one-dimensional res- 
onant photonic crystals. In order to achieve a complete understanding of relationships 
between optical and geometrical characteristics of these structures, one needs a general 
analytical method, which could be applied to a variety of RPCs independently of a partic- 
ular spatial dependence of their dielectric function. In our previous paper, Ref.y, such a 
method was developed for studying reflection, transmission, and absorption spectra of one- 
dimensional RFC structures. The present paper presents a general analytical approach to 
calculating the band structure of these materials. This method allows one to analyze disper- 
sion laws of one-dimensional RFCs with an arbitrary form of the periodic modulation of the 
dielectric function, propagation angle of the electromagnetic wave and its polarization state. 
This method naturally incorporates resonant excitations of the medium into the theory. For 
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concreteness, we assume that the resonances are related to exciton states in QWs, but this 
approach can be easily applied to any other type of the resonance excitations. In Sections ITTl 
and, particulary, in subsection 111 Bl we show that by choosing an appropriate basis for rep- 
resenting transfer matrices one can significantly simplify the analysis of the band structure 
and dispersion laws of electromagnetic excitations in the structures under consideration. In 
Sections IV and V we illustrate the power of the approach by applying it to a number of 
examples, some of which deal with already well studied situations (passive one-dimensional 
photonic crystals, optical lattices, MQWs with a piece-wise spatial dependence of the dielec- 
tric function), while the others demonstrate the possibilities of our approach in obtaining 
results that could not have been obtained by other methods. 

II. THE (a, /)-REPRESENTATION OF THE TRANSFER MATRIX 
A. Transfer matrix approach in resonant photonic crystals 

A propagation of the electromagnetic wave in structures under discussion is governed by 
the Maxwell equation 



where z axis is chosen along the growth direction and Pexc is the excitonic contribution to 
the polarization, which can be presented in the following form: 



where $^(-2) = $(-2 — is the envelope wave function of an exciton localized in the m-th 
QW. The summation in Eq. (j2)) is performed over all QWs and z^ are the positions of their 
centers. We assume that the distance between the consecutive wells, d = Zm+i — Zm, coincides 
with the period of the spatial modulation of the dielectric function, e{z + d) = e{z). Also, we 
assume that the profile of the dielectric function is symmetrical with respect to the position 
of the center of the QW, e{zm + z) = e{zm — z) (see Fig. H]). We restrict ourselves to the 
consideration of Is states of heavy-hole excitons and neglect their in-plane dispersion. The 
dipole moment of these excitons lies in the plane of the well, therefore, they can only interact 
with the in-plane component of the electric field, Ej^. This fact is reflected in Eq. (j2I), where 
only this component of the field is included. 




(1) 
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FIG. 1: An example of the modulation of the dielectric function, e(z) = e + Ae{z) (solid line). The 
dashed vertical lines show the positions of the centers of QWs. 

The frequency dependence of the excitonic response is given by the function, xi^)^ defined 



as 



a 



(3) 



ujq — uj — t'y 

where Uq is the exciton resonance frequency, 7 is the non-radiative decay rate of the exciton, 
and a is the microscopic exciton-light coupling parameter, proportional to the dipole moment 
of the electron - heavy hole transition. 

Eq. can be analyzed by presenting the direction of the electric field as a vector sum 
of two mutually perpendicular polarizations. These, so called s- and p-polarizations, define 
two possible eigen directions, and accordingly, the respective electric fields satisfy two inde- 
pendent equations, which can be considered separately. In the structures modulated in the 
2;-direction, normal modes of the system are characterized by the wave vector perpendicular 
to the modulation (growth) direction, which is a conserving quantity. The direction of this 
vector can be also used to define the eigen polarizations of the waves. 

The absence of an overlap of the exciton wave functions localized in different QWs makes 
it possible to use the transfer matrix technique for solving Maxwell equations, Eq. (Q). The 
general idea of this technique is to characterize the electric field of the wave with a certain 
polarization by a two-component column vector, \c{z)) whose change along the structure is 
governed by a 2 x 2 matrix. The main property of this matrix can be described as following. 
If \c{zi)), \c{z2)), and \c{z3)) characterize the field at points Zi, Z2, and 2:3, respectively, and 
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if T{zi, Zj) satisfies the equation \c{zi)) = T{zi, Zj)\c{zj)), then T{zi, z^) = T{zi, Z2)T{z2, z^,). 



The relations between the values of the field and the vector \c{z)) as well as the form of the 
matrix T{zi,Zj) depend on the choice of the particular representation of the field. Several 
possible representations will be discussed in the subsequent sections of the paper for both 
polarizations. 

1. S -polarization 

In the s-polarized wave the electric field E is perpendicular to the plane formed by the 
direction of z axis and in-plane wave vector, k. Accordingly, it can be presented in the form 



where p is the coordinate in the (x, ?/)-plane and = ek x is a unit polarization vector, 
defined in terms of unit vectors in z- and k-directions. For E[z) we obtain an integro- 
differential equation 



where k1{z) = uj'^e{z)/c^ — k'^. 

We first consider the transfer matrix corresponding to the propagation of the field across 
one elementary cell of the structure, i.e. from point + at the left boundary of the cell, 
inside it, to the point z+ + just outside of the right boundary of the cell. The complete 
transfer matrix connecting the field at the right boundary of the entire system with the field 
at its left boundary can then be obtained as a product of the transfer matrices for each 
cell. Considering only one cell, we can, without any loss of generality, choose the origin of 
our coordinate system coinciding with a position of the QW in the cell under consideration. 
When z coordinate is confined to a single cell, the summation over QWs and the index 
enumerating them in Eq. (0) can be dropped. 

The term with the exciton polarization at the r.h.s. of Eq. (0) can be considered as an 
inhomogeneity in a second order differential equation 



E{z,p) =esE{z)e'' 



(4) 





d^E{z) 
dz^ 



+ kI{z)E = T{z) 



(6) 
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whose general solution can be written in the forn>sa 



E{z) =cih{z)+C2h2{z) + {Gi.J^){z) 



where 



{Gi<J^){z) 



dz'J-'{z' 



Mz')h2{z)-h,{z)h2{z') 



(7) 



(8) 



W{h^,h2]Z' 

Here we have introduced hi^2{,z), a pair of linearly independent solutions of the homogeneous 
equation 



(fE{z) 



+ Ki(z)E = 0, 



(9) 



and W{hi, h2', z) = hih'2 — h[h2 is the Wronskian of these solutions. For the case under 
consideration the Wronskian does not depend on z and will be denoted by Wh in what 
follows. 

Alternatively, a solution of Eq. ^ can be presented as 



E{z) = ci{z)hi{z) + C2{z)h2{z), 
E'{z) = c^{z)h[iz)+C2iz)h'2iz). 



(10) 



In the regions outside of the QW, where ^{z) = 0, functions c{z) remain constants, making 
representations given by Eq. ((Zj) and Eq. (fTO|) equivalent, but q(z) change when z traverses 
from one boundary of the QW to the other. Our goal now is to use Eq. (0) along with 
Eq. (jS)) in order to relate the values of Ci{z) at the right boundary of the QW to their values 
at the left boundary. Using these equations we can obtain the expression for the electric 
field for points to the right from the QW as 



E{z) =hi 
+h2 



ci + X 2 — (^1*^1 + "^aV^s 

C2 - X 2 (^iV^l + 



(11) 



where ^\^2 are "projections" of the solutions /ii^2 onto the exciton states 



1 



d^^{z!)hx;2{z) 



JQW 

and the modified excitonic response function x can be presented as 

X 



(12) 



X 



(13) 
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where 

Auj = a I dz<!>{z){Gi.<^>){z) (14) 
Jqw 

gives the radiative shift of the exciton frequency in the photonic crystal. Eq. (fT^ is a 
generahzation of the well-known expression for the radiative shift in MQW structures with 
a homogeneous dielectric function .^^i^^i^^ 

Comparing Eq. (fTTj) with Eq. (fTUj) . and taking into account that the coefficients q do not 
change outside of the QW, we can find a relation between the coefficients Cj determined at 
the two inner boundaries of the elementary cell, z = + 0, and z = z^ — 




with matrix Th given as 



(15) 



T„^l + *-^h^' I, (16) 



where 1 is the unit matrix. It should be noted that Eq. (fTB|) is valid for an arbitrary form 
of the exciton envelope wave function and spatial modulation of the dielectric function. 

Eq. can be significantly simplified if we use freedom in the choice of the functions 
hi{z). In the case when e{z) is invariant with respect to the mirror refiection relatively to the 
center of the QW, hi^2 can be chosen to have a definite parity, with one of them being even 
with respect to the center of the QW, and the other being odd. With such a choice of these 
functions we can always turn one of v?i,2 to zeroA- It should be understood, however, that hi 
and h2 are not actual normal modes of the photonic crystal. The latter must be defined as 
solutions of the appropriate boundary problem, and do not have to be even or odd functions 
with respect to the center of the elementary cell. Moreover, due to the Bloch theorem the 
modes of a photonic crystal can have a definite parity only at specific frequencies that are 
naturally identified with boundaries of the forbidden gap in the spectrum. However, all 
normal modes of the photonic crystal can be represented as superpositions of the even and 
odd solutions hi^2- A formal discussion of the relation between the functions hi and /12 and 
the normal modes can be found in Ref. |27| where the similar approach has been used for the 
analysis of spectral properties of a Schrodinger equation with a periodic potential. 

Obviously, fixing the parity of the solutions does not determine functions hi and /i2 
uniquely since a multiplication by a constant determined by initial conditions does not 



change their symmetry. It can be shown, however, that results obtained for observable 
quantities, such as band structure, do not depend on this ambiguity. We will present our 
general results in a form independent on the choice of initial conditions. However, for 
discussions of particular examples it is more convenient to fix them in the form 

hiO) = l, h[iO) = 0, 

(17) 

/i2(0) = 0, h'^iO) = 1, 

for which Wh = 1. 

Let us assume that hi is the even solution, then = ^ and the expression for the matrix 
Th simphfies to 

Tft = 1 + S.5. , (18) 

where Qs = k,s{z+) is the value of k,s{z) at the boundary of the elementary cell and 



s. = -n^f-^. (19) 



Substitution of x yields 

Ss = V— (20) 

uj — ujq — Auj + «7 



where Tg is the radiative decay rate. 



(21) 



The function Ss{uj)-, which we shall refer to as exciton susceptibility, describes the contribu- 
tion of the exciton-light interaction to the optical properties of the structures under consid- 
eration. For instance, the resonant absorption of light occurs at the frequency u}q = ujq + Au; 
determined by the pole of ^^(ti;). We will treat ujq as an experimentally accessible exciton 
frequency, which along with the radiative decay rate Fg can be measured in optical exper- 
iments with a single QW. In what follows we will drop the tilde from ujq and assume that 
the radiative shift is included in this parameter. 

Coefficients Cj in Eq. (|15|) are the two components of the vector |c) characterizing the 
field in the transfer matrix formalism while matrix given by Eq. is the transfer 

matrix itself. The functions hi{z), which we used to derive Eqs. (|15|) and (jl8|) form a basis 
for this particular representation of the transfer matrix. Using this basis we are able to 
obtain a simple and convenient expression for the transfer matrix describing the evolution 
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of the electric field inside a single elementary cell. However, in order to obtain the transfer 
matrix through the period of the structure, we have to relate the field and its derivative at 
points to the left {z^ — 0) and to the right (z+ + 0) of the boundary between two elementary 
cells. Since there is no discontinuity of the dielectric constant at this point, it might appear 
that all the continuity conditions would be satisfied automatically. One should remember, 
however, that functions hi{z) are different in different elementary cells. For two adjacent 
m-th and (m + l)-th cells, the respective functions, h^{z) and h'^^^{z), are related to each 
other as h^(z) = h^^^^z — d), and in order to establish connection between vectors |c™) and 
\c"^^^) at the boundary between these cells, one would need to express h'^~^'^{z) in the basis 
of h^{z). This problem can be solved, but it turns out to be more convenient to convert our 
transfer matrices to a more conventional basis of plane waves 

E{z) = E+e'"^' + E.e-'"^' (22) 

using the conversion rule, Eq. ()A9|) . derived in the Appendix. The resulting transfer matrix 
describing the evolution of the field across an entire elementary cell can then be presented 
in the form 

\(a]-ra)l2 af ) 
which, as will be seen shortly is particularly convenient for analysis of the spectrum of the 
system under consideration. This representation is extensively used in the present paper 
and, in what follows, we will refer to it as (a, /)-representation. The parameters of this 
representation, a and /, are defined as 

a = 92, f = 9i- iSs92, , ^ 

(24) 

a = ^2) i = 9\^^Ss9l, 

where 




2. p -'polarization 

A representation similar to Eq. ()23p for a transfer matrix describing p-polarized light can 
also be obtained along the similar lines. An important difference is that the electric field in 
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the p-polarized waves is not transverse (V ■ E 7^ 0) and, therefore, satisfies a more complex 
system of two differential equations. For conventional photonic crystals the p-polarized 
waves are conveniently described in terms of the magnetic field (see e.g., Ref. 1^. However, 
the presence of the dipole active excitations significantly reduces the convenience of this 
approach because in order to close the equation for this field, one would have to express the 
interaction with excitons in terms of the magnetic field. While it can be done, we find it 
more convenient to continue working with electric field, E. 

The electric field for the p-polarization can be presented in the following form: 

E{z, p) = [e^,E,{z) + e-,E-,{z)] e'^f . (26) 

Amplitudes Ex^z{z) satisfy the system of equations 

d'^Ea: dEz 2 



dz'^ dz 



ik— h k(z)E^ 



-X(^)5^$™(2;) J dz<^rniz')E^{z'), (27) 



dE^ 
dz 

where Kp{z) = uj'^e{z)/(?. Deriving these equations we again have explicitly taken into 
account that only the in-plane component of the electric field interacts with the heavy-hole 
excitons. Solving the second equation with respect to Ez we obtain the closed equation for 

EJz) 



d 
dz 



\ .dE, 
p{z 



dz 



-X{uj)^^2^<^^{z) / dz<^^{z')E^{z'), 

m 

where p{z) = k^{z) /[k^{z) — k"^]. This function is related to the local angle of propagation 
of the wave, 9{z): p{z) = 1/ cos6'(2;). 

The derivation of the transfer matrix for the p-polarized field follows exactly the same 
steps as in the previous subsection with one important change. The Wronskian of two 
solutions of the "homogeneous" (that is with % = 0) version of Eq. (PHj) depends on the 
coordinate z^ This, however, does not present a major problem, because if one takes into 
account that this dependence can be presented in the form 

W{z) = W,^, (29) 
12 



one can immediately reproduce all the main results of the previous subsection. In particular, 
the transfer matrix in the basis of plane waves can again be written down in the (a, /)- 
representation, Eq. (|23|). The parameters of this representation are given by the same 
Eqs. (I24|) and (j25j) . where, however, one has to consider functions hi^2 as even and odd 
solutions of the "homogeneous" version of Eq. (|28|). Also, the function Sp, which retains the 
same form as in Eq. ()20|). is now characterized by a modified radiative shift of the exciton 
frequency and the decay rate. The expression for the frequency shift is obtained by replacing 
the expression given by Eq. (jH)) with 

= Wfl{z+) [ ^^'^^^'^ [hi{z)h2{z) - hr{z)h2{z')] (30) 

in Eq. p4|l for the radiative shift. The rate of exciton's decay into the p-polalized light is 
given by 

r, = ^^. (31) 

where qp = Kp{z+). One can notice that these expressions and corresponding expressions for 
the s-polarized wave coincide in the case of normal propagation, i.e. when k = 0. 



B. General properties of the (a, /)-representation 

The fact that the transfer matrices for both polarizations allow for the (a, /)- 
representation is not a coincidence but is related to a well-known feature of the Maxwell 
equations. To demonstrate it let us consider Eq. in a particular case corresponding to 
the passive structure (x = 0). If E{z) is a solution to this equation, its complex conjugate 
E*{z) is also a solution. Having this in mind one can show that that 

^ [p{z) {EE*' - E'E*)] = 0. (32) 

Now, representing the electric field in the form of Eq. (j22|l we can find that 

^^[qpPiz){\E^\'-\E_\')]=0. (33) 

This relation expresses conservation of the flux of the Poynting vector through a plane 
perpendicular to the 2;-axis and represents the fact that there are no sources or drains of the 
energy in the system. The similar expression for the s-polarized wave can be obtained by 
substitution qp qg, p{z) 1. It follows from Eq. that the transfer matrices through 
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the period for both s and p polarizations must preserve the combination — l-E'-P, 

which means that the transfer matrix written in the basis of plane waves belongs to SU (1,1) 
group.— A general form of an element of this group is^ 



(34) 



where T12 are complex numbers and |Tip — = 1. Matrices allowing for the (a, /)- 
representation correspond to a particular case of SU{1, 1) matrices, for which T2 is purely 
imaginary. It can be shown that such matrices describe structures that posses a mirror 
symmetry. Indeed, a transfer matrix of such a structure satisfies the following condition 
(JxTa^ = T~^, where is the Pauli matrix. Substitution of Eq. ()34|) into this equation gives 
T2* = — T2. The converse statement can also be verified. This property of transfer matrices 
can also be proven in the general case of matrices corresponding to Maxwell equations with 
external polarization (0) and (j28|) . but the proof is rather technical and we do not provide 
it here. 

Owing to the (a, /) -represent at ion, the description of structures with mirror symmetry 
can be significantly simplified. At the same time, they demonstrate a rich variety of inter- 
esting phenomena, and, therefore these structures attract a great deal of attention (see e.g. 
Ref. I29). If such a structure is built of blocks that have mirror symmetry by themselves then 
the transfer matrix through the entire structure can be written in terms of matrix elements 
describing the individual blocks. Indeed, let us consider a structure with the period BAB 
where the transfer matrices of the blocks A and B have the form T(ai, /i) and T(a2, /2), 
respectively, in the (a, /)-representation. Then it can be shown that the transfer matrix 
through the period is 

T(a2, /2)T(ai, /i)T(a2, /2) = r(a, /), (35) 

where 

a = 0102/2 - ^ (02/2 - 02/2) , 

/i - ^^^^ 

/ = /2/ia2 + — (02/2 - 02/2) • 

We will illustrate the "multiplication rule", Eq. (j35p . considering a simple but very im- 
portant for the rest of the paper example of an elementary cell in an optical lattice. In this 
system, the block A corresponds to an active element, for instance, a QW, and two blocks 
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B describe barriers, which have the same refractive index as block A. The blocks B are 
described by a transfer matrix given by Tb{(f)b) = diag [exp(i0b), exp(— z^;,)]. In this case one 
has 02 = /2 = exp(i0f,/2) in Eq. and thus 

%{<Pb)na, f)TM = nae^f"', fe^'^"). (37) 

If, however, one needs to take into account the difference between refractive indices of 
the blocks B, and the block A, a new type of transfer matrix products emerges. Introduc- 
ing matrix Tp, describing propagation of light across a boundary between two media with 
different indexes of refraction 

where p is the Fresnel reflection coefficient^ (see below Eq. ()49p). we can describe propa- 
gation of light through BA and AB interfaces and across the block A using the following 
combination of the transfer matrices: Tp{p)~^T{a, f)Tp{p). An important property of the 
(a, /)-representation is that this product can be expressed in the following form: 

Tpip)-'T{a, f)Tp{p) = Y^Tia + ap, f - fp). (39) 

The real factor in Eq. (IHUj) can be incorporated into a and / due to the useful relation (with 
real A) 

XT{a,f) = T{Xa,f) = T{a,Xf). (40) 

Eqs. (jHHjl . (jSZ|), and will be extensively used throughout the paper for obtaining 
the (a, /)-representation of various transfer matrices. This method is often more practical 
than solving corresponding differential equations. Thereby, it is interesting to note that the 
parameters of the (a, /)-representation are essentially boundary values of solutions of the 
corresponding Cauchy problem. 

III. THE DISPERSION EQUATION IN RESONANT PHOTONIC CRYSTALS: 
THE METHOD OF ANALYSIS AND THE BAND STRUCTURE 

A. Dispersion equation in the (a, /)-representation 

The dispersion equation characterizing normal modes (polaritons) in an infinite periodic 
one-dimensional structure can be expressed in terms of elements of a transfer-matrix across 

15 



one period of the structure^ 

cosfsTrf = ^TrT, (41) 

where K is the Bloch wave number, d is the period of the structure, and the transfer matrix is 
assumed to be written in the plane wave basis. If this structure possesses a mirror symmetry, 
one can show that the condition det T = 1 can be presented as fd + fa = 2. Using this 
identity we can rewrite the polariton dispersion law in the form 

^) = 5J(a)5R(/), (42) 

where 

sft(a) = {a + a)/2. (43) 
Equivalently, this dispersion equation can be presented as 

sin^ (^) = %aMf), (44) 

where 

5>(a) = (a - a)/2i. (45) 

If the exciton susceptibility S" is a real function (no homogeneous broadening of excitons), 
the a coincides with a conjugated, d = a*. In this case, 9ft(a) and '^{a) are equivalent to 
regular real or imaginary parts of a respectively: 3?(a) = Re(a) and Q{a) = Im(a). When 
homogeneous broadening is taken into account, Eq. remains valid, but identification of 
3?(a) and 53(a) with Re(a) and Im(a) can no longer be made. In this case, one has to use 
definitions given in Eqs. (j43|l and (j45|l . 

If the homogeneous broadening can be neglected, the r.h.s of Eq. (jl^ and Eq. (jl^ are 
real valued. Therefore, we can analyze the band structure in the vicinity of the exciton 
resonance using the notion that for allowed bands (real K) these expressions should be 
positive and less than unity, while for the band gaps (complex K) they should be negative 
or greater than unity. Thus, there are two types of conditions determining the boundary of 
the bands. From Eq. ()42|1 we have that the band boundary occurs when either 9ft(a)9ft(/) 
or 3?(a)3?(/) — 1 changes sign, and Eq. defines as the boundaries those frequencies at 
which the change of sign occurs for expressions Q{a)'^{f) or Q{a)Q{f) — 1. One can show 
that these two pairs of conditions are equivalent to each other. More precisely, expression 
3fJ(a)3fJ(/) changes sign at the same frequencies as expression Q{a)Q{f) — 1. The same is 
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true for the pair 3fJ(a)3?(/) — 1 and ^5(a)3(/). Thus, all the band boundaries can be found 
by considering changes of sign of expressions 3fJ(a)3ft(/) and '^{a)Q{f), with negative values 
corresponding to the band gaps. The factorized form of right-hand sides of the polariton 
dispersion equation, Eqs. (j42j) and (j44|) . drastically simplifies the analysis of the spectrum, 
as it will be seen in the subsequent sub-sections. 

The factorization of the dispersion equations, Eqs. ()42|) and fj44|) becomes possible only 
owing to the (a, /)-representation of the transfer-matrix, which makes this representation 
particularly suitable for studying the band structure and dispersion laws of polaritons in 
resonant photonic crystals. At the same time, this representation is not very convenient 
for calculating relfection/transmission spectra of finite structures because of a cumbersome 
relation between the parameters of a single layer transfer matrix and a matrix describing 
the entire structure. This problem, however, can be conveniently solved with the help of 



23. It is useful to establish 



a different representation introduced in our recent paper, Ref. 
a direct relationship between these two representations. First we notice that an arbitrary 
matrix of the form can be written in the form of a transfer matrix describing propa- 
gation of a wave across a single QW [see Eq. below]. This can be done by introducing 
an effective optical width of one period of the structure, (p, and an effective excitonic sus- 
ceptibility, S. The relation between the parameters of the (a, /)-representation and these 
effective parameters can be found by comparing Eqs. and 

S=^{af-af), e^^ = -. (46) 

The usefulness of this relation is twofold. First, it shows that for any resonant photonic 
crystal the transfer matrix through the period can be represented as a transfer matrix 
through a conventional MQW structure. Second, Eq. (j46|) can be used for a consideration 
of reflection and transmission spectra of a structure described by a transfer matrix in (a, /)- 
representation. This is where the introduction of S and and then, consecutively, the 



representation the transfer matrix in terms of parameters 6 and (3 of Ref. |23| becomes very 
convenient. However, this representation misses the possibility to factorize the dispersion 
equation, so it is not very suitable for a constructive analysis of the dispersion law. One can 
see, therefore, that the representations of the transfer matrix in terms of parameters (a, /) 
and S and compliment each other in a sense that each of them is suited best for its own 
set of problems. 
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FIG. 2: The periodic structure built of two blocks. Vertical dashed lines show the boundary of the 
elementary cell. The angles of propagation inside the blocks are related by Snell's law. 



B. MQW structures with a simple elementary cell 

As it was mentioned in Introduction, MQW structures present one example of RPCs, 
in which effects due to reflection of light from well-barrier interfaces coexist with effects 
due to the resonant light-exciton coupling. The spatial profile of the refractive index 
in these structures is particulary simple, and is described by a piece-wise constant func- 
tion [Fig.©]. Various properties of these structures have been studied in a number of 
publication o^'^1^'^1^^1^^1^^1^^1^'^ , including papers Refs|20||2l||23| dealing particularly with disper- 
sion laws and the band structure of their normal modes. Nevertheless, we find it useful to 
consider this case in details in this paper because it gives a clear illustration of using the 
(a, /) -represent at ion for the analysis of complex dispersion equations. 

We start by considering two particular cases: a passive photonic crystal (no resonances) 
and an optical lattice (no refractive index mismatch). This will give us useful benchmarks for 
discussing the general situation. In the latter case we will neglect the exciton homogeneous 
broadening, which allows us to identify operations 3? and with taking regular real and 
imaginary parts of the respective expressions. The same is obviously true for the case of 
passive structures. 



1. Passive photonic crystals 

Let us consider a structure built of a periodic sequence of two blocks characterized by 
the widths dh^^ and the indices of refraction nf,^^- To emphasize the mirror symmetry we 
choose the elementary cell as shown in Fig. |21 The transfer matrix through the period of 
the structure has the form 

T = Tl'''T;^T^TpTl'\ (47) 
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where 

e"*"^" 

(pb^w = ^nh^wdb,w cos 9b,w/c and 9b,w are the angles of propagation of the wave. 



(48) 



The scattering of the electromagnetic wave at the interface between different blocks de- 
pends on both the angle of incidence of the wave and its polarization state. These effects 
are described by Fresnel coefficients ps and Pp, 

riw cos 6yj — rib cos 6b 



P 



s 



riw cos 6yj + Hb cos 9b ' 



cos 6b - nb cos 6^ . 

Pp = y^^) 

riw cos 6b + Ub cos 6^ 

for s and p polarizations, respectively. Below we denote the Fresnel coefficients simply by 
p having in mind that for a particular polarization only one of these expressions should be 
used. 

Using Eqs. (|!?7j) and we can find parameters of the (a, /)-representation of the transfer 
matrix, which in this case take the following form 

« = " 2 ^ / = e^"^" - Pe^"^" ' (50) 

1 — P^ 

where t± = {(pb ± 0w)/2ci;. Now we can use general dispersion equations, Eq. (jl^ and 
fl44|) to analyze the structure of the allowed and forbidden bands of this system. As we 
discussed in Section Fill Al there are two types of forbidden bands, determined by conditions 
Re (a)Re (/) < 0, and, Im (a)Im (/) < respectively. If r+ and r_ are incommensurate, 
these band gaps alternate along the frequency axis and can be classified as either odd (first, 
third, etc.) or even (second, forth, etc.). 

Since the coefficients of the dispersion equation in the passive photonic crystals do not 
contain any singularities, its right hand side can only change the sign by passing through zero. 
Thus, the band boundaries are situated at the frequencies where either Re (a) Re (/) = or 
Im (a)Im (/) = 0. For instance, if p > 0, the boundaries of odd band gaps are determined 
by equations 

Re [a(fi+)] = 0, Re [fifl^)] = 0, (51) 

where f2_ and i7+ correspond to the lower and higher frequency boundaries of a given band 
gap respectively. The explicit form of these equations can be obtained using definitions. 
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Eq. dSni): 

cos(r2_ri ) — pcos(^7_r_) = 0, 

(52) 

cos(i7+r+) + pcos(i7+r_) = 0. 
In the simplest case, when the layers have the same optical width, one has r_ = and the 
positions of the edges of the forbidden gap are Ur{l ± 2 arcsin(p)/7r), where 

UJrT+ = ^. (53) 

While this case gives a convenient reference point, however, having in mind applications to 
MQW structures, an opposite situation, when the optical widths of the layers are different, 
is of more interest. Generally, as one can see from Eq. the boundaries of the gap are 
situated asymmetrically with respect to ujr. Under the assumption of narrow gaps, which is 
fulfilled for a not very strong contrast of the refraction indexes, and for angles not too close 
to the angle of the total internal reflection, the positions of the boundaries are given by 

n^ = u:Jl±-^^^^^). (54) 



7r(l ± pcos(j)w) 



2. Semiconductor optical lattice 



In this case we assume that all layers in the structure have the same index of refraction, 
n, but there are dipole active excitations in QWs. The propagation of light through a QW 
is described by the transfer matrix of the form 

fe^'/'^n-iS) -iS \ 
Tw=\ , 55 

y iS e~'^-{l + iS)j 

where 0^ = und^, cos 9u,/c. The parameters of the (a, /)-representation of are 

a = g*0™/2^ f ^ ei<A»/2(l _ ^gy (56) 

The excitonic contribution to the scattering of light is described by 

S = (57) 

The radiative decay rate, Fq, depends on the angle of incidence. As follows from Eqs. (pT|l 
and (ISH), for structures with homogeneous dielectric function these dependencies for different 
polarizations ar o^^i^^i^^ 

r^^) = To/ cos r;,^^ = To cos (58) 
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The transfer matrix through the period of the structure is 

T = tI''^T^tI''^ . (59) 

Taking into account the multiphcation rule (jHTj) . we can find parameters of the {a, f)- 
representation for the entire matrix ()59|) and the respective dispersion equatior ^'^'^^ 



cos^ (^^~) ~ cosc<jr_|_(cosc<jr+ + Ssmur^). (60) 

Unhke the case of the passive photonic crystal, the coefficients of the (a, /)-representation 
now contain a singularity at the frequency of the exciton resonance (we neglect the homo- 
geneous broadening of the excitons here). Therefore, the r.h.s. of the dispersion equation 
(jHn|l can change its sign by passing not only through zero, but also through infinity. The 
latter happens at cuq, which becomes one of the boundaries of the band gap associated with 
the exciton resonance. In general case, the width of this and all other allowed and forbidden 
bands are proportional to F, which is many orders of magnitude smaller than ujq. Therefore, 
usually, exciton related modifications of the photon dispersion law in the optical lattice are 
negligible. This situation changes, however, if we require that this singularity is canceled by 
the term cos(c<jr+), which happens if cos(ct;oT+) = 0. This is equivalent to the condition for 
the Bragg resonance, when the half-wavelength of the radiation at the exciton frequency is 
equal to an odd multiplier of the period of the structure 

uJoT+ = - + vrn. (61) 

The spectrum of such structures is characterized by a much wider band gap with the width 
equal to 

Ar = 2./^ = 2.0^°^ (62) 
and is indicative of enhanced coupling between light and QW excitons. 

3. MQ W structures with a mismatch of the indices of refraction 

Now we are ready to discuss the general case of MQW structures with the contrast of 
the refractive indices. The transfer matrix for this structure is obtained from Eq. by 
taking into account the scattering at the interfaces between the QWs and the barriers, 

T = Tl'^T;^T^TpTl'\ (63) 
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The dispersion equation following from the (a, /) -represent at ion of this transfer matrix has 
the form 

cos^ (^^^ = -l-^Re{apc)[Reifpc) + ^lm(apc)], (64) 

where apc and fpc are calculated for a passive multilayer structure and are given by 
Eqs. 

Similarly to the case of the optical lattice, the structure of the spectrum near the exciton 
frequency is complicated by the singular character of the excitonic susceptibility (see Fig. OJ, 
which gives rise to new branches of exciton related collective excitations^^ and respective 
band gaps. We, however, will again focus on Bragg structures, characterized by the can- 
cellation of the excitonic singularity. As a result, an anomalously broad band gap in the 
vicinity of the exciton frequency is formed. As one can see from Eq. (jl^ . such a cancelation 
occurs when either Re(apc(t^o)) = or Im(apc(co'o)) = 0. Both these cases imply that the 
exciton frequency coincides with a boundary of one of the photonic band gaps. Most exper- 
iments with these systems tend to deal with structures having as short period as possible. 
Therefore, we only consider the case, when uq = where is the boundary of the lowest 
band gap of the respective photonic crystal, determined by the first of equations (j^T|) . The 
exciton related band gap in this case appears at the frequencies where the r.h.s. of Eq. fl64|) 
is negative. 

Assuming a smallness of the gap we can expand apc and fpc near the frequencies fl± 
Re[apc(c^)] = {n+ - iu)t+, Re[/pc(cu)] = ((]_ - iu)t_, (65) 
and write the equation for the boundaries of the forbidden gaps in the form 

{to - n+) (iu-n.)- To^^^^^^ = 0, (66) 



where 



^Re[apc{u)] 



(67) 



In Eqs. fj65|) we explicitly have taken into account the negative sign of these derivatives. In 
Eq. (jUBj) the radiative decay rate, Tq, is defined in Eqs. fISHjl should be taken according to 
the polarization of the wave and the angle of propagation. 

Using these approximations we can present the boundaries of the polariton band gap as 

u;± = Qc± ^A, (68) 
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FIG. 3: The figure plots the r.h.s. of Eq. (|64() scaled by 10^ for a structure slightly detuned from 
the Bragg resonance as a function of frequency. The material parameters are chosen to be close 
to typical parameters of GaAs/AlajGai-^^As structures: p = 0.03, Fq = 60 /ieV, ujq = 1.5 eV. The 
frequencies where this expression is negative or exceeds unity (not shown in this plot) correspond 
to band gaps. The vertical line shows the position of the exciton frequency. 

where VL^ = + fi+)/2 is the center of the gap and A is its width. The expression for A 
can be written in the following form: 



A 



A%c + A2 



(69) 



and is equal to "Pythagorean sum" of the widths of the passive photonic gap, Apc 
— and the modified excitonic gap. 



A^ 



To 



4:lm{apc) 



(70) 



t_ ' 1 - P 

The last formula in this equation is obtained by using approximate expressions lm{apc) ~ 
1 + p and t_ ^ r+(l — p) which are valid provided that the optical width of the wells is very 
small compared to the width of the barriers. 
The equation 

coo = n+ (71) 

generalizes the condition for the Bragg resonance given by Eq. for optical lattices. 
Indeed, normal modes in photonic crystals are characterized by the Bloch wave number Kpc, 
rather than by the wave number of a homogeneous medium un/c. For the band boundaries 



23 



uo = Q± one has Kpc{^±)d = vr, so that ujq = is equivalent to Kpc{uJo)d = vr, which is 
a direct generahzation of Eq. ()61|) . expressed in terms of the appropriate wave number. 

While the results presented have been obtained using the assumption that p > (that 
is, for the normal propagation, > n^) they remain valid in the opposite case due to the 
symmetry of the transfer matrix under the transformations p —p and ctpc ^ fpc- Iii 
other words, when we have the opposite relation between and Ub all the arguments used 
above can be repeated with the mirror reflection of the frequency axis with respect to the 
center of the photonic gap. In particular, the Bragg resonance occurs when the exciton 
frequency coincides with the left (low frequency) edge of the photonic band gap. 

4. Angular dependence of the band structure 

In previous publications on the band structure of MQW systems^! only modes propagat- 
ing along the growth direction of the structure were considered. In this paper, thanks to 
a general nature of our approach, we can consider waves propagating at an arbitrary angle 
treating both s- and p- polarizations of the electromagnetic waves on equal footing. The 
foundation for this consideration is laid by the results of the previous subsections of the 
paper, where the expressions for the parameters defining the band structure were derived in 
terms of Fresnel coefficients, Eqs. fl49j) . The latter contains all the information about angular 
and polarization dependencies of resonant as well as refractive index contrast contributions 
to the band structure. 

It is clear from Eq. that the Bragg frequency defined by the generalized Bragg 
condition, Eq. (f7T|) . depends on the propagation angle of the wave. This dependence is 
presented in Fig. |3] for both s- and p-polarized modes. Due to the narrowness of the QWs 
the position of the Bragg resonance approximately follows the renormalization of the optical 
width of the period of the structure oc cos{9b) for both polarizations. This effect can be used 
to tune the position of the band gap of the structure. Indeed, if one changes the position 
of the exciton frequency, by, for instance, using quantum confined Stark effect,— the system 
can be tuned back to the Bragg resonance by changing the propagation angle of the wave. 
Assuming that the structure remains tuned to the Bragg condition with the changing angle 
we can consider the angular dependence of the width of the polariton band gap. One can see 
from Fig. |3J where this width is presented by vertical bars, that its angular dependence is 
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FIG. 4: Dependence of the Bragg resonance frequency lobI^t [see Eq. ()53() ] on the angle of prop- 
agation 0^ measured in the barriers. The material parameters are the same as in Fig. |31 Bold and 
thin lines correspond to p— and s-polarizations, respectively. The vertical lines (the error bars) 
show the forbidden gap for each polarization. For better visibility the gap is scaled by the factor 
of 5. 

different for different polarizations. For the s-polarization the width monotonously increases 
when the angle increases, while for the p-polarization it, first, decreases, reaches its minimum 
at the Brewster's angle determined by the equation sin6'b = n^j \/nl + nl^ (it corresponds 
to 6b ^ 0.8 in Fig. |3]), where the only contribution to the gap is due to the exciton-light 
interaction, and then starts increasing. It should be noted that at the Brewster's angle the 
Fresnel coefficient pp changes its sign. Thus, as it has been pointed above, the condition of the 
Bragg resonance for the p-polarized wave propagating at angles larger than the Brewster's 
angle is formulated as equality of the exciton frequency and the lower frequency edge of the 
photonic band gap for the p-polarized wave. 

Structures with the fixed exciton frequency and the period can be tuned to the Bragg 
resonance at a single angle only. For waves propagating at other angles the Bragg condition 
is violated, and the structure becomes detuned from the resonance. In the case of small 
detunings, it is natural to refer to such systems as off-Bragg systems. Thus, the discussion 
of the angular dependence of the band structure, naturally involves considering properties 
of the off-Bragg structures. 

When the exciton frequency is not tuned to one has to take into account the sin- 
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gularity of the exciton susceptibility when looking for the boundaries between allowed and 
forbidden bands. This singularity results in a narrow band gap situated between zeros of 
Eq. (jH)), i.e. where siv?{Kd/2) < 0. This band gap occupies the region between ujq and 
ujQ + VLi, where 

= \r+5Kl (72) 

and 5 = ujq — VL^ ^ 0. For slightly off-Bragg structures, the width of this addition to the 
band gap is small, and we will neglect its existence in the future discussions. Thus, the 
band boundaries are determined by zeros of r.h.s. of Eq. (jl^ . and as the result the band 
structure is determined by four frequencies: uq and the roots of the equation 



iuj — Vt, 



0. (73) 



The band structure is characterized by two band gaps and one allowed band between 
them. If we put these four frequencies in ascending order, the band gaps will lie between 
the first and second pairs of frequencies separated by a transparency window (see Fig. Oj). 
The exact order of the band boundaries depends on relation between uq and If uq < 
then the band gaps are between and ujq and between f2+ and uj'_^, where 

S . 1 



UJ. 



± 



± J(Apc-Sy + Al, (74) 



while the allowed band is between c^o and Thus, detuning the exciton resonance fre- 
quency away from leads to the appearance of the transparency window between ujq and 
fi+ inside the band gap obtained for the Bragg case Uq = and to a slight modification of 
the external edges of the gap. This result is in a qualitative agreement with the earlier anal- 
ysis of the off-resonant MQW structures.— The transparency window would manifest itself 
in optical spectra as a dip in the reflection coefficient. A similar dip was actually observed in 

urn 

Refs. 191 )231)391 where structures with the period satisfying the non-modified Bragg condition, 
Eq. ()UT|) . were considered. Since the well and barrier materials had different, albeit close, 
values of the refractive indexes, the real Bragg period should have been determined from 
Eq. ()71|) . and the structures considered in those papers were actually slightly off-Bragg. It 
should be mentioned, however, that the detuning from the real Bragg condition could be 
not the only cause for the observed dip in reflection. The inhomogeneous broadening of the 
excitons could also contribute to this effect.-^ 
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FIG. 5: The dependence of the band gap structure of the s-polarized wave on the angle of propa- 
gation 9h measured in the barriers. The material parameters are the same as in Fig. |31 The dashed 
regions correspond to the forbidden gaps. The structure is assumed to be tuned to the Bragg 
resonance at normal propagation. 

Now we can describe how the band structure of our system evolves with the changing 
propagation angle. Let us assume that the system is at the Bragg resonance for the normal 
incidence. Increasing the angle, we detune it away from the resonance. As a result, two band 
gaps, one adjacent to the exciton frequency and the other one detached from it, separated by 
the transparency window, appear. Fig.^lshows the angular dependence of the band structure 
for the case of the s-polarization. An interesting phenomenon is seen to occur at an angle 
at which = uj'_^_. In this case, the width of the band gap detached from the exciton 
frequency turns to zero (band gap collapse). This collapse is a specific feature of RPCs 
and is absent in both purely passive structures and optical lattices. An important property 
of the gap adjacent to uq is its weak polarization dependence. This can be understood if 
one notices that any difference between polarizations could only appear at angles where 
ub — Ar, Apc". The width of the adjacent gap for these angles can be found from 

Eq. dZH). It is equal to 

With the accuracy up to the terms ~ Aadj Apc/ {ub — ujq), it is the same for both polariza- 
tions. The property of the omnidirectional reflectivity, i.e. a resonant reflection at all angles 
of incidence, ensues from these results. It follows from the notion that a usual scattering 
problem is set up for a structure embedded in a medium with essentially lower index of 
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refraction (air or vacuum). As the result, the angles of propagation inside the structure can- 
not exceed the angle of the internal reflection 9c at the boundary between the structure and 
surrounding medium. Therefore, for all angles of incidence the range of frequencies between 
ujq and uj'_{9() corresponds to the forbidden gap and, therefore, to a resonant reflection. The 
similar effect has been considered for the case of passive photonic crystals in a number of 
publications.'^i'^^iiSi^ Here we want to emphasize the feature specific for the RPC If the an- 
gle of total internal reflection is small (for the air-GaAs interface dc ~ 0.28), we can describe 
the change of the edges of the photonic band gap by a simple renormalization of the optical 
width of the period, VL± Q±/ cos 6*. It results in the width of the region corresponding to 
the omnidirectional reflection in the form 

1 A2 

A ■ ^ - =^ (76) 

""""^ 2 fi+sin2^,-2Apc' ^ ^ 

where we have assumed that the photonic forbidden gap is not too wide, i.e. f2+ sin^ 9^ > 
2Apc- One can see that the presence of QWs essentially weakens the condition of the 
omnidirectionality in comparison with passive photonic crystals. 

C. General modulation of the dielectric function 

1. Resonant photonic crystals 

The formalism developed in this paper allows us to analyze the band structure in 
RPCs with an arbitrary periodic modulation of the dielectric function. Using the (a, /)- 
representation for transfer matrices and expressions for the parameters a and / given by 
Eqs. (im), we can derive the dispersion equation for such a structure in the form similar to 
Eq. (jUl): 

or Eq. dm) 

2 



sm 



(f) = IKi^.) + S,/4(.H-)| , (78) 

These dispersion equations can be used to describe waves of all polarizations by choosing 
appropriate for a given polarization parameters S, q, and hi 2. The band structure of a 
passive photonic crystal can be obtained from these equations by setting x = 0. The band 
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FIG. 6: Dependence oilin(Kd) (solid lines, left axes) and wave functions and its derivative (dashed 
and dotted lines, respectively, right axes) on the frequency for a passive structure {S = 0). The 
non-zero values of lm{Kd) correspond to the band gaps. The index of refraction is n{z) = 3 + 
cos^'^(7rz/2). (a) Kd is the solution of Eq. ((77|) . /ii(a'), the dotted line depicts h2{uj). (b) Kd has 
been found from Eq. ((75|) . the dashed line represents h2{uj), the dotted line depicts h[{io)/5. One 
can see how zeroes of the wave function and its derivative determine consecutive band boundaries. 



boundaries, Q±, in this case are determined by equations 



/ii(2;+,fi_) = 0, /i^(2;+,fi+) = 0, 



(79) 



and 



0, 



(80) 



where in the argument denote that the functions hi^2 are solutions of corresponding 
homogenous equations [Eqs. (0) or with x = 0] for = Qzp respectively. Thus the 
zeroes of even and the maxima of odd solutions given by Eqs. (f?^ and (jHUJ determine 
the boundaries of the consecutive band gaps in the band structure of the system under 
consideration (see Fig. IHI). 

The analysis of the full dispersion equations, with the resonance terms restored, repeats 
the analysis of the previous subsection, where one needs to make substitutions apc = g2 
and fpc = 91- In particular, the parameters t± of the equation for the boundaries of the 
forbidden gap, Eq. are defined in terms of boundary values of /i,2 as 



dh'2in^ 



(81) 
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The general structure of the spectrum is similar to the one discussed in the case of the 
piece-wise modulation of the refractive index, and we can generalize the condition for the 
Bragg resonance, when the exciton related allowed band collapses, and the spectrum of 
our system in the vicinity of the exciton frequency consists of a single wide band gap. We 
again require that the singularity of the excitonic susceptibility at Uq is canceled by the 
first term /i2(-2+) in Eq. (f77j) . This corresponds to the condition ujq = where f2+ is the 
high frequency boundary of the respective photonic gap. From the dispersion equation ()79p 
it follows again that the Bragg condition can be casted in the form Kpc{ujQ)d = vr. One 
can rewrite the last equation in yet another form emphasizing the role of the phase in the 
formation of Bragg structures. Let us consider the effective optical width of the period of 
the structure, [see Eq. (j46|) ]. At the point where /ig = one finds that a/a = —1 yielding 
4>{^l+) = TT. Thus, the phase form of the Bragg condition, which is particularly convenient 
for practical calculations, becomes (f){uJo) = tt. 

We can also generalize the expression for the width of the band gap in the case of the 
Bragg resonance. Performing the expansion similar to Eqs. we can obtain the equation 
for boundaries of the forbidden gap: 

(u - n+) I iu - n_ - ^^^^ I = 0, (82) 

y UJ -UjQ J 

with Ap given by the first part of Eq. ()70|). The analysis of Eq. ()82j) completely repeats that 
of Eq. (jUBj) . In particular, the width of the forbidden gap is determined by the Pythagorean 
sum of the photonic and excitonic contributions, 

= Ale + ^r, (83) 

where Apc = — and expression for the excitonic contribution, Ar, can be presented 
as 

f|£y.^I£!* and f^V = ££!M(±t> (84) 

for s- and p-polarizations, respectively. Here we have introduced 

«S = ^ / dze{z)E^^f, (85) 

where Ei|p is the s- and p-polarized electric fields corresponding to the even mode of the 
photonic crystal. This expression represents the energy associated with the even mode of 
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the photonic crystal concentrated in a single elementary cell. The r.h.p. of Eqs. ()84|) can be 
shown to be proportional to u'qw/u''^\ where u~lw is the energy of the even mode concentrated 
within the QW. This means that the contribution of the exciton-light interaction into band 
structure of the resonant photonic crystal depends on the distribution of the energy of the 
electric field inside the elementary cell at the frequency corresponding to The greater 
amount of energy stored within the well the larger the role of the excitonic effects is. 

2. Off-Bmgg structures 

Our analysis so far was limited to consideration of the band boundaries and the properties 
of the band gaps in the resonant photonic crystals. To complete the consideration of the 
normal modes in these structures we need to discuss solutions of Eqs. (f?^ and (f?^ in the 
allowed bands of the spectrum. These solutions determine dispersion laws uj{K) of the 
photonic crystal polaritons. In the case of the Bragg structures and for frequencies in the 
vicinity of the band edges, the dispersion laws can be presented in the following form: 

u±{K) =n,± ^VA2 + (2(^-7^5)2, (86) 

where Kb = vr/rf is the Bloch wave number corresponding to the boundary of the first 
Brillouin zone and C, = d/t^t_. These two branches with positive and negative masses, 
m± = ±2h'^A/('^, lie in the bands that are above and below the forbidden gap, respectively. 
In addition to these branches there is one more dispersionless branch at the exciton frequency 

UJb{K) = Uq. 

For the systems detuned from the exact Bragg resonance, the positions of the band 
boundaries and properties of the band gaps can be determined by repeating the arguments 
from the previous sub-section; we do not reproduce it here. We however, complement that 
analysis by considering modifications of the dispersion laws of the polariton branches for the 
off-Bragg structures. If the detuning is not too large, \6\ = \uJo — ^ Ar, the dispersion 
laws of the original (present in the Bragg structure) upper and lower polariton branches 
near the band edges can be written in the following form 

uj±{K) =uj'^ + !f{K~ Ki,f/2m±, (87) 

where u'^ are given by Eq. ([7^1 with parameters defined by generalized equations (fTHjl . (jH21), 
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and ()84|) . The renormalized mass parameters m± are defined as 



m±{6) = m± I 1 



A2 



± 
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) 



(88) 



Comparing this result with Eq. ()86|) one can notice that detuning from the Bragg resonance 
results in two main modifications of the dispersion laws: first, the position of the band 
boundaries have changed and, second, the magnitudes of the masses of the upper and the 
lower polariton branches are no longer equal to each other. 

The most dramatic changes occur, however, with the third, originally dispersionless 
branch. In off-Bragg structures this branch acquires dispersion, which, of course, agrees 
with the opening of the allowed band in the vicinity of the exciton frequency. The disper- 
sion law characterizing this band is given by 



This branch corresponds to excitations dubbed "Braggaritons" in Ref. The mass of this 
mode, given by = A'^/2(6, is very sensitive to the amount of detuning from the Bragg 
condition and can be effectively controlled, for instance, by the electric field via the quantum 
confined Stark effect. This property of the " Braggariton" branch invited proposal to use it 
for slowing, stopping and storing light in BMQW structures.— 

3. Homogenous broadening 

We finish our consideration of dispersion properties of RPCs with an arbitrary periodic 
distribution of the refractive index by discussing effects due to the exciton homogeneous 
broadening. In the presence of dissipation the concept of band gaps becomes ill-defined 
because the imaginary part of the Bloch wave-number, generally speaking, is not zero at all 
frequencies. Nevertheless, the properties of this imaginary part are of great interest, since 
they determine the spectral and transport properties of the structures under consideration in 
the vicinity of what would have been band boundaries in the absence of dissipation. In order 
to access these properties we have to consider solutions of Eq. (f77j) . Taking into account 
that in the spectral region of interest, the real part of K is close to n/d, we look for solutions 
to this equation in the form 




T 



(89) 



Kd = IT + iX. 



(90) 
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FIG. 7: Dependence of lm{Kd) on frequency for different structures in the vicinity of the first for- 
bidden gap when the exciton homogeneous broadening is taken into account. Dashed line represents 
the passive structure from Fig. El Solid line shows a Bragg MQW structure with a homogeneous 
dielectric function. Dotted line represents the structure with combination of QWs and the smooth 
modulation of the dielectric function. The characteristic feature is a divergence of the penetration 
length (ImK)~^ at the exciton frequency. 

The unknown parameter A here has a simple physical meaning: the inverse of its real part 
determines a characteristic length on which the amplitude of the mode decays. The same 
length determines the penetration depth of incident radiation, therefore it is often called a 
penetration length. Assuming that A is small we can find an approximate expression for it 
in the form: 



An important feature described by this expression is that A vanishes at the frequency cor- 
responding to the right edge of the photonic band gap, signifying the divergence of the 
respective penetration length. The similar effect takes place in MQW^° structures without 
the refractive index contrast. The important difference is that in the latter case the diver- 
gence occurs at the exciton frequency that lies at the center of the forbidden gap rather 
than at the band's boundary. Fig. [7| shows the comparison of the solutions of the disper- 
sion equation for different structures with the exciton homogenous broadening taken into 
account. 




(91) 
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D. MQW structures with complex elementary cells 



In this section we illustrate the application of the developed technique to structures with 
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In that paper the mismatch 



complex elementary cells, which were first studied in Ref. 
between refractive indexes of different elements of the structure had to be neglected because 
the combination of the complex elementary cell and the spatial modulation of the refrac- 
tive index turned out to be an insurmountable obstacle for the standard transfer-matrix 
approach. In this section we show that the approach developed in the present paper al- 
lows us to overcome the technical difficulties associated with the consideration of a complex 
structure, and to generalize the results of Ref. I2] for the case of structures with modulated 
refractive index. 

We will focus here on one particular example, when the elementary cell consists of two 
QWs with different exciton frequencies, and UJ2, located half a period apart from each 
other (see Fig. IH]). It was found in Ref. |2l| that despite the presence of two different exciton 
frequencies, it is still possible to generalize the notion of the Bragg resonance for such systems 
and design structures whose spectrum would consist of only two polariton branches separated 
by a wide band gap. However, unlike the case of structures with a simple elementary cell, 
the formation of such a wide band gap requires not only the period of the structure to 
have a certain value, but also imposes a condition on the spectral separation between the 
excitonic frequencies of the wells constituting the elementary cell. When both generalized 
Bragg conditions are fulfilled for such a structure, the width of the polariton band gap, Acs, 
becomes larger than that in the case of structures with a simple elementary cell: 

Acs = V2Ar. (92) 

This broadening of the gap by almost 40% refiects a possibility to effectively strengthen the 
exciton-photon interaction by increasing the density of QWs in the structure. 

In this section we consider how the mismatch of the refractive indices affects the spectral 
properties of such structures. We simplify our consideration assuming that all elements of 
the structure except the wells with the exciton frequency uj2 have the same refractive index. 
Formally, the dispersion equation describing modes propagating in such a structure in the 



normal direction has been obtained in Ref. 
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by more conventional methods, but that 



equation turned out to be too cumbersome to allow for any non-numerical analysis. Using 
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FIG. 8: The periodic structure with two QWs (dark rectangles) in the elementary cell. Dash 
lines show the boundaries of the elementary cell having the mirror symmetry. The QW with the 
exciton frequency uj2 is assumed to have the index of refraction different from other elements of 
the structure. 

(a, /) -represent at ion of the transfer matrices, we show here that the dispersion equation can 
be rewritten in a much more transparent form with a factorized right hand side. 

To apply the technique developed in the present paper it is necessary to choose the 
elementary cell with the explicit mirror symmetry. It can be done as shown in Fig. |H| The 
problem of determining the transfer matrices through the right and left halves of the QW 
can be resolved in the following way. We note that the QW transfer matrix determined 
by Eq. can be written in a factorized form 

= T6(0^/2)T;„Tf,(0^/2), (93) 

where Ty^ is derived from the expression for by setting = 0. After this factorization 
the transfer matrix through the elementary cell takes the form similar to Eq. (j35p 



T = \/fi'^ n{<p, + 0^/2) r-iri,2) TpT,{<p, + 0^/2) ^/fi'K (94) 

The difference between the indices of refraction is again taken into account by introducing 
the special transfer matrix Tp, which describes propagation of the wave across the interface 
and contains respective Fresnel coefficients. The transfer matrices through different QWs, 
Tw , are obtained from Eq. (j55j) by substituting different excitonic susceptibilities 

-51,2 = (95) 



The matrix square root \/ T^^ can be found to be equal to 



1 - iSi/2 -iSi/2 



Ti'^= I ' I . (96) 

tSi/2 1 + iSi/2 
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Now, we can apply Eqs. (j^ . (jH7|) and to derive (a, /)-representation for T. The 
dispersion equation following from this representation has a relatively simple form: 

1 = [Re(a) + Si Im(a)] [Re(6) + ^2 Im(6)] , (97) 

/ 1 - P 

where a = exp(z0+) — pexp(z0_), 6 = exp(20+) + pexp(z0_) [compare with Eqs. and 
0± = oj{2dbn + (i^n ± du]n2)/2c. 

As usual, the positions of band boundaries are determined by frequencies at which the 
r.h.s. of Eq. (jHTjl changes its sign, which can occur either by passing this expression through 
zero or through infinity. The other set of boundaries, which we do not consider here can be 
obtained as before by converting Eq. (jU7|) to the form with sin^ (-^j on the left. We can show 
that in the immediate proximity of the exciton frequencies, there are six band boundaries, 
which define three band gaps separated by two allowed bands. We again will be interested 
in finding conditions, under which these allowed bands collapse resulting in formation of one 
continuous band gap. As we already know from the previously considered examples, the 
collapse of the allowed bands associated with the singularities of the exciton susceptibility, 
S", takes place when the exciton frequency coincides with one of the band boundaries of a 
certain effective structure. In the case of a simple elementary cell this effective structure was 
a respective passive photonic crystal, however, in the case under consideration, the effective 
structure is comprised of one of the sublattices of our system. More precisely, in order to 
cancel the exciton singularity at, for example, cui, we have to require that uji coincides with 
the zero of the expression inside the second parenthesis on the r.h.s. of Eq. (jU7|) . This 
zero, obviously gives the position of some band boundaries of a structure made up of the 
sub-lattice 2. Similarly, in order to cancel the singularity at uj2 we need to make it equal to 
one of the band boundaries originating from the sub-lattice 1, particulary that which turns 
to zero the expression in the first parenthesis of the r.h.s of Eq. (j^Tj) . Let us assume for 
concreteness that uj2> ooi-, then more specifically we have to require that uJi is equal to the 
low-frequency boundary of the band gap associated with polariton branches in sub-lattice 2, 
and UJ2 should coincide with the high-frequency boundary of the band gap originating from 
polariton branches in sublattice 1. Let us denote frequencies of these two band boundaries as 
ojf2i where the lower index shows from which sublattice the respective boundary originates, 
and ± denotes high- and low-frequency boundaries, respectively. The frequencies 10^2 can 
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be found to be equal to 



(98) 



where fl^p are determined by Eqs. (j52j) . and for simplicity, we assumed the same value 
of the effective exciton-light interaction for both sublattices, i.e. we neglected a possible 
difference between exciton radiative decay rates in different wells. We also neglected the 
renormalization of Fq due to the spatial modulation of the dielectric function [see Eqs. fl70p 
and ()84|) ]. The generalization of these expressions is straightforward but does not lead to 
anything new while making the resulting expressions much more cumbersome. 

Let us, first, consider the case when the index of refraction of the wells of the second type 
is higher, i.e. p > 0. The conditions of absence of transparency windows inside the gap (the 
Bragg resonance conditions) are in this case formulated as = uj2 and u;^ = oji. Resolving 
these equations with respect to the exciton frequencies one obtains 



which generalize the results of Ref. to the case of structures with the refractive index 
mismatch. Substituting these results into the expression for the width of the band gap, 
which is equal to A = cu^ — cuf , we obtain a familiar result 



The case of p < is described by the same formulas with substitution Apc —Apc, which 
changes the Bragg condition [the second equation of Eqs. (jSHl)] but leaves the width of the 
gap, Eq. (|l(J(jp . the same. 

Thus, taking into account the mismatch of the indices of refraction in the system with 
a complex elementary cell leads to modification of the Bragg condition in comparison with 
what one has in the case without the mismatch. This modification has a different character 
than in structures with a simple cell. The mean value of the exciton frequencies stays at 
the center of the forbidden gap and only the difference between the frequencies must be 
corrected to take into account the mismatch. The resultant gap is, similar to the case of 
simple cell, again the Pythagorean sum of the photonic and the excitonic contributions. 



2 




(99) 




(100) 
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IV. CONCLUSION 



We considered the band structure associated with normal excitations of one-dimensional 
resonant photonic crystals using as a particular example the exciton polaritons multiple 
quantum well structures. Properties of normal modes of such structures are determined 
by the interplay of two mechanisms: the interaction of light with internal resonances of 
constituting materials, such as excitons in MQW case, and multiple reflection of light due 
to periodic spatial modulation of the index of refraction of the structure. The necessity to 
take into account both mechanisms of the light-matter interaction in these structures on 
equal footing makes the analytical description of their properties more complicated than in 
the case of regular passive photonic crystals. In this paper we developed a powerful method 
of analyzing band structure and dispersive properties of such structures, which allowed us 
to obtain analytical solutions for problems that could not be solved by other approaches. 

This method is based on two key features. First, we construct transfer matrices us- 
ing solutions of the appropriate Cauchy (rather than boundary value) problems as a basis. 
This makes calculations of the transfer matrices more convenient since it allows for avoid- 
ing difficult questions about eigenmodes of the system under consideration. Second, we 
use a special representation of the transfer matrix valid for structures with an elementary 
cell possessing a mirror symmetry. This representation automatically yields the dispersion 
equation of normal modes in a factorized form, which drastically simplifies the analysis of 
the spectrum. In particular, using this representation we were able to analyze dispersion 
properties of certain type of resonant photonic crystals with complex elementary cells. We 
also believe that besides situations explicitly considered in the paper, the suggested method 
allows for effective treatment of other problems such as exciton luminescence from MQW 
photonic crystals, or for analysis of more complicated systems, e.g. quantum graphs. ^'''^^ An 
important benefit of the developed approach consists also in the fact that polarization and 
angular dependencies of the elements of transfer matrices appear in this approach only via 
respective Fresnel coefficients. This enabled us to obtain band structures of modes of both 
s- and p-polarizations in MQW based photonic crystals from essentially the same dispersion 
equation and easily analyze their angular dependencies. 

We illustrated the developed method by applying it to a number of problems, some of 
which were considered previously in the framework of other approaches, while others were 
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analyzed in this paper for the first time. In particular, we considered polariton spectrum 
in several systems: MQW structures with the mismatch of the indices of refraction of the 
barriers and wells materials, MQW with an arbitrary periodic modulation of the dielectric 
function, and MQW with a complex elementary cell consisting of two equidistant QWs 
characterized by different exciton frequencies. We showed that in all these cases one can 
construct a structure, in which multiple photon and polariton bands can be collapsed in 
just two polariton branches separated by an anomalously broad stop band. Such a band 
structure is characteristic of so called Bragg structures, which in the case of MQWs were 
studied in a number of previous publications. We showed here that in the case of RPCs with 
a simple elementary cell the condition for formation of such a spectrum can be formulated as 
a requirement for the exciton frequency to coincide with one of the edges of the respective 
passive photonic crystal. Formulated in this form this condition demonstrates a direct 
connection with earlier formulations of the Bragg resonance in MQW structures without the 
refractive index contrast. In the case of crystals with complex elementary cells, we showed 
that our method allows for taking into account the mismatch of the refractive indices between 
different elements of the structure and found the generalized Bragg conditions for this system 
as well. In the last case, this condition consists of two equations: one relating the period of 
the structure to the exciton frequencies of the constituent elements and the other specifying 
a relation between those frequencies. 

One of the surprising results is that for all these structures the width of the polariton band 
gap in the case of the Bragg resonance can be presented in the form of a Pythagorean sum of 
the photonic and the excitonic contributions, = Ap,-. + Ap. The excitonic contribution 
in the presence of the refractive index modulation, however, is modified compared to the 
case of simple optical lattices. It is interesting that this modification can be presented in the 
form of a ratio of the electric field energy concentrated inside QW to the the total energy 
stored in the entire elementary cell of the structure. 
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APPENDIX A: TRANSFER MATRIX FORMALISM 



The transfer matrix technique has been reviewed in a number of pubhcations. Here, 
however, we use this technique in a somewhat unusual form. Therefore, we find it useful to 
recall some important results. While the general idea of the transfer matrix establishing a 
linear connection between values of a function representing a solution of a linear differential 
equation at two different values of the argument is inherent to all formulations, the explicit 
form of the transfer matrix and its properties depend upon a choice of the representation 
of the matrix. As a basic representation one often chooses a 2 x 2 matrix operating on a 
two-dimensional vector, whose components are the value of a function, which represents a 
solution at a given point, and its derivative. The formal reason for this choice is the fact 
that a linear equation of a second order can always be naturally reformulated as a system 
of two first order equations for the function and its derivative: 



Another approach to constructing transfer matrices is more convenient for problems in- 
volving scattering of waves. This problem is most naturally described by specifying an 
amplitude of an incident wave and relating to it amplitudes of transmitted and reflected 
waves. In this approach, the field is sought in the form e***^ + re~^'^^ at the one side of the 
system, and in the form te^'^^ at the the other side. More generally, the field is represented 
in the form 



and the transfer matrix relates the amplitudes at the left and the right boundaries 




(Al) 



E{z) = a+e*''^ + a_e" 
E{z + d) = a'+e'"' + a'_e 



(A2) 




(A3) 



As follows from Eqs. ()A2|) . the transfer matrices written in these two bases are related by a 
similarity transformation, T = W~^{z)T^{z + d, z)We{z), where 




(A4) 
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Another important way to introduce the transfer matrix (see e.g. Ref. 1491 ) is based on 
a possibihty of a representation of the solution of the Maxwell equation as a sum of two 
linearly independent functions, /i,2(-2), with modulated amplitudes 



E{z) = ci{z)hi{z) + C2{z)h2{z), 
E'{z) = c,{z)h[{z) + C2{z)h'2{z). 



(A5) 



Now, the transfer matrix gives the relation between the amplitudes at different points, 
usually a period of a structure apart. 



Ci{z + d) 
C2{z + d) 



Th{d) 



ci{z) 



C2(Z 



(A6) 



From Eq. ()A5|) it follows that 



nid) = W^\z + d)T^{z + d, z)Wh{z), 



(A7) 



where Wh{z) is the Wronsky matrix 



Wh{z) 



hi{z) h2{z) 
h[{z) h',{z) 



(A8) 



Eq. ()A8|) allows one to derive a relation between transfer matrices obtained for a different 
choice of the basis functions hi^2- 

The relation between the transfer matrices written in the bases of plane waves and a pair 
of linearly independent functions can be written as 



T = M{d/2)Th{d)M-\-d/2), 
where M{z) = W-\0)Wh{z). Using Eqs. (jMl), IM), and (jMl) one obtaines 



(A9) 



( 



hi{z) + 
hAz)- 



h\{z) 



h2{z) + 



h',{z)\ 



IK ^ IK ^ 



(A 10) 



IK IK / 

All types of transfer matrices have their advantages and disadvantages. For instance, a 
basis of a pair of linearly independent functions can be most naturally related to solutions of 
the Cauchy problem for the respective differential equation. Such a basis is easy to find but 
it is inconvenient for the solution of scattering problems. The fact that, generally, it is not 
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related to by a similarity transformation [see Eq. ()A7|) ] makes it also difficult to consider 
the spectral problems using this basis. Indeed, the knowledge of the transfer matrix itself 
is not sufficient for solving this problem, and one has to separately consider a boundary 
value problem for amplitudes Ci^2- On the contrary, the basis of plane waves is the most 
suitable for solving this particular problem since the transfer matrix itself contains all needed 
information. In this basis, the difficulty of having to solve the boundary value problem is 
moved to the process of finding the transfer matrix itself. Therefore, the conversion rule 
from one basis to another presented in Eq. ()A9|) is very useful for practical calculations as 
it is demonstrated in our paper. 
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